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Abstract 

Increased  control  demands  in  applications  including  high  speed  milling  and  hybrid  motor  design  have  led  to 
the  utilization  of  magnetostrictive  transducers  operating  in  hysteretic  and  nonlinear  regimes.  To  achieve  the  high 
performance  capabilities  of  these  transducers,  models  and  control  laws  must  accommodate  the  nonlinear  dynamics 
in  a  manner  which  is  robust  with  regard  to  system  inputs  and  facilitates  real-time  implementation.  This  necessitates 
the  development  of  models  and  control  algorithms  which  utilize  known  physics  to  the  degree  possible,  are  low-order, 
and  are  easily  updated  to  accommodate  changing  operating  conditions  such  as  temperature.  We  consider  here  the 
development  of  nonlinear  adaptive  identification  techniques  for  low-order,  energy-based  hysteresis  models  having 
nonlinear  parameterizations.  We  illustrate  the  techniques  in  the  context  of  magnetostrictive  transducers  but  they 
are  sufficiently  general  to  be  employed  for  a  number  of  commonly  used  smart  materials  including  piezoceramics, 
magnetostrictives  and  shape  memory  alloys.  The  performance  of  the  resulting  nonlinear  identification  algorithms  are 
illustrated  through  numerical  examples. 

Index  Terms:  Magnetostrictive  materials,  hysteresis,  constitutive  nonlinearities,  nonlinear  parameterization, 
adaptive  estimation 


1  Introduction 

Piezoceramic,  magnetostrictive  and  shape  memory  alloy  materials  are  being  increasingly  considered  for  use  in  high 
performance  applications  due  to  their  unique  transduction  capabilities.  Piezoceramic  transducers  are  lightweight, 
exhibit  both  actuator  and  sensor  capabilities,  are  highly  effective  for  broadband  and  high  frequency  transduction,  and 
have  extremely  high  set  point  accuracy.  These  properties  are  utilized  in  applications  ranging  from  the  design  of  PZT- 
based  stages  for  nanopositioning  in  an  atomic  force  microscope  (AFM)  [11,  24]  to  inertial  sensors  in  accelerometers 
[10].  Magnetostrictive  transducers  also  exhibit  broadband  transduction  capabilities,  actuator  and  sensor  capabilities, 
and  the  ability  to  generate  large  forces  while  maintaining  highly  accurate  set  point  placement.  This  has  led  to  their 
use  as  actuators  in  industrial  applications  such  as  high  speed  milling  and  highly  sensitive  sensors  for  applications 
such  as  torque  sensing  in  a  steering  system  [6,  21].  Shape  memory  alloys  have  the  largest  work  output  densities  of  the 
three  classes  of  materials  but  are  inherently  slower  than  piezoceramics  and  magnetostrictive  transducers  due  to  the 
cooling  required  as  a  prelude  to  heat-induced  phase  transformations.  The  large  strain  capabilities  of  SMA  actuators 
has  motivated  their  use  in  applications  such  as  active  tendon  design  to  attenuate  earthquake  and  wind-induced 
vibrations  in  civil  structures  [25],  microgripper  design  for  biomedical  instruments  [9],  and  large  force  generation  for 
reconfigurable  airfoil  design  [13]. 

While  these  materials  provide  unique  transduction  capabilities,  the  advantageous  performance  attributes  come 
at  the  cost  of  hysteresis  and  constitutive  nonlinearities  throughout  most  of  the  drive  range.  Furthermore,  these 
inherent  material  properties  are  dependent  upon  environmental  conditions  and  system  inputs  such  as  temperature, 
drive  frequency,  and  input  or  residual  stresses  -  see  [16,  23]  for  quantification  of  these  dependencies  in  the  magne¬ 
tostrictive  transducer  compounds  Terfenol  and  Galfenol.  To  achieve  the  transducer  capabilities  in  moderate  to  high 
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Figure  1:  (a)  Prototypical  magnetostrictive  transducer  employed  for  high  speed  high  accuracy  milling,  (b)  Hysteretic 
field-magnetization  relation  as  reported  in  [7]. 


drive  regimes,  it  is  necessary  to  quantify  the  hysteresis  and  constitutive  nonlinearities  in  a  manner  which  is  suffi¬ 
ciently  accurate  to  characterize  the  fundamental  nonlinear  and  hysteretic  material  behavior  and  sufficiently  efficient 
to  facilitate  transducer  design  and  real-time  implementation.  A  fundamental  step  in  the  model  construction  entails 
the  estimation  of  model  parameters  given  various  system  inputs.  Furthermore,  it  is  often  necessary  to  adaptively 
update  parameters  to  accommodate  the  inherent  changes  in  material  behavior  which  occur  as  environmental  con¬ 
ditions  change.  For  example,  ohmic  heating  in  a  magnetostrictive  transducer  can  significantly  affect  the  hysteretic 
field-displacement  relations  exhibited  by  the  transducer  whereas  Joule  heating  in  an  electrostrictive  transducer  can 
eliminate  the  hysteresis  if  temperature  changes  are  sufficiently  large.  While  these  environmental  changes  often  occur 
on  a  relatively  slow  time  scale,  their  effects  on  the  model  must  be  accommodated  to  maintain  the  fecundity  of  the 
transducer.  In  this  paper,  we  develop  nonlinear  parameter  estimation  techniques  in  the  context  of  energy-based  mod¬ 
els  which  have  been  employed  to  quantify  the  hysteresis  and  constitutive  nonlinearities  inherent  to  magnetostrictive, 
piezoceramic,  and  shape  memory  alloy  transducers. 

To  focus  the  discussion,  we  consider  a  prototypical  magnetostrictive  transducer  employed  for  applications  including 
high  speed,  high  accuracy  milling  as  depicted  in  Figure  la.  Input  forces  and  displacements  are  provides  by  the 
Terfenol-D  rod  in  response  to  fields  generated  by  the  surrounding  solenoid  due  to  an  applied  current.  As  detailed 
in  [4,  7,  8],  these  transducer  designs  have  the  capacity  for  generating  broadband,  high  force  responses.  However, 
they  also  exhibit  varying  degrees  of  hysteresis  and  constitutive  nonlinearities  in  the  relations  between  input  fields 
H  and  the  generated  magnetization  M  and  strains  e  as  illustrated  in  Figure  lb  as  detailed  in  [7].  Moreover,  these 
nonlinear  and  hysteretic  relations  are  temperature-dependent  and  thus  change  as  the  transducer  heats  during  the 
milling  process.  This  necessitates  the  development  of  adaptive  parameter  estimation  techniques  to  update  models  to 
maintain  the  specified  cutting  accuracy  during  the  milling  process. 

While  this  application  is  quite  specific,  it  illustrates  a  number  of  the  issues  shared  by  magnetostrictive,  piezo¬ 
ceramic  and  shape  memory  alloy  transducers  operating  at  moderate  to  high  drive  levels  in  variable  temperature, 
variable  stress  operating  conditions.  Furthermore,  it  will  be  demonstrated  that  the  modeling  framework  and  nonlin¬ 
ear  parameter  estimation  techniques  are  sufficiently  general  to  permit  direct  extension  to  a  wide  range  of  piezoceramic 
and  SMA  devices. 

There  exist  a  number  of  techniques  for  modeling  hysteresis  in  ferromagnetic  and  magnetostrictive  materials 
including  micromagnetic  relations  [3,  12],  Preisach  representations  [1,  22,  32],  domain  wall  models  [4,  7,  8,  14,  15], 
and  homogenized  free  energy  relations  [27].  Whereas  micromagnetic  theory  addresses  fundamental  mechanisms 
associated  with  magnetostriction,  its  complexity  presently  precludes  its  use  in  real-time  control  systems.  Preisach 
models  have  the  advantage  of  a  general  and  rigorous  mathematical  framework  and  can  be  linearly  parameterized  in 
terms  of  coefficients  to  be  identified,  thus  permitting  the  use  of  linear  adaptive  estimation  and  control  algorithms  [32] . 
These  advantages  are  offset,  however,  by  the  large  number  of  nonphysical  parameters  required  for  general  operating 
regimes  and  the  theoretical  extensions  required  to  accommodate  frequency,  stress  and  temperature  dependence.  The 
domain  wall  models  are  low-order  and  have  a  small  number  of  parameters  (e.g.,  5-8)  which  can  be  correlated  with 
physical  properties  of  the  data.  Moreover,  the  methodology  is  sufficiently  general  to  encompass  ferroelectric  and 
piezoceramic  materials  [28,  30]  and  ferroelastic  and  shape  memory  compounds  [19]  so  it  provides  a  unified  framework 
for  quantifying  hysteresis  and  constitutive  nonlinearities  in  a  broad  range  of  ferroic  compounds  [29].  For  these  reasons, 
we  employ  a  domain  wall  model  to  quantify  the  field-strain  behavior  depicted  in  Figure  lb.  However,  the  models 
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exhibit  a  nonlinear  dependence  on  parameters  which  motivates  the  development  of  nonlinear  adaptive  estimation 
techniques  to  accommodate  the  temperature-dependent  hysteresis  behavior  exhibited  by  the  materials. 

In  this  paper,  we  extend  the  techniques  of  [17,  18]  to  obtain  nonlinear  adaptive  estimation  laws  for  updating 
parameters  in  the  domain  wall  models.  Among  the  advantages  of  the  parameter  estimation  algorithms  presented 
here  are  the  properties  that  they  can  be  implemented  during  the  operation  of  the  actuator  or  sensor  and  require  no 
additional  data  collection  as  required  for  offline  least  squares  algorithms.  One  possibility  is  to  employ  the  parameter 
estimation  algorithm  during  actuator  operation  to  identify  parameters  during  the  initialization  of  the  milling  process 
as  a  prelude  to  control.  Since  the  behavior  of  the  materials  is  affected  by  changing  operating  conditions,  such  as 
temperature,  the  estimation  algorithm  can  be  used  to  re-identify  parameters  in  the  transducer  model  to  reflect 
changes  in  the  material  behavior. 

A  significant  difficulty  in  developing  nonlinear  parameter  adaptation  laws  is  the  fact  that  gradient  update  methods 
are  not  always  sufficient  for  estimating  nonlinearly  occurring  parameters.  To  illustrate,  consider  an  error  model  of 
the  form 

e  =  -Ae  +  /0, 9)  -  /(</>,  9) 

where  /  is  a  nonlinear,  scalar-valued  function,  </>  denotes  measured  states,  9  and  9  respectively  denote  a  nonlinearly 
occurring  vector  of  parameters  and  their  estimate,  A  >  0  is  a  scalar,  and  e  is  the  error  between  desired  and  measure 
trajectories.  Consider  the  gradient  update  law 

9  =  eVfd. 

With  the  standard  Lyapunov  function,  V  =  \{e2  +  92),  where  9  =  9  —  9,  we  see  that 


V  =  -Ae2 


/(<M)-/(<M)  +  0V/g 


(1) 


If  e  <  0,  it  is  necessary  that  V/g(0  —  9)  <  9)  —  f(cf>,  9)  which  is  satisfied  if  /  is  convex.  If  e  >  0,  then  V  <  0  is 

ensured  if  /  is  concave.  Hence,  the  gradient  method  does  not  ensure  stability  for  all  9.  Moreover,  a  gradient  method 
applied  to  a  nonlinear  parameterized  system  may  not  only  be  insufficient  but  may  lead  to  instability.  The  methods 
we  discuss  here  do  not  strictly  rely  on  a  gradient  rule  but  differ  depending  on  the  sign  of  the  error  as  motivated 

by  (!)■ 

In  Section  2,  we  summarize  the  hysteresis  and  transducer  models  used  to  characterize  the  magnetostrictive 
transducer  depicted  in  Figure  la.  In  Section  3  we  review  the  nonlinear  adaptive  method  for  the  scalar  case  given 
in  [17,  18]  and  present  the  extension  of  this  theory  to  a  vector  system  as  is  required  for  the  Terfenol-D  transducer. 
Section  4  provides  numerical  examples  of  a  nonlinear  parameter  estimation  for  both  the  scalar  and  vector  cases. 


2  Actuator  Model 

In  this  section,  we  construct  models  for  the  magnetostrictive  transducer  which  can  be  employed  in  the  nonlinear 
parameter  estimation  algorithm  given  in  Section  3.  We  develop  these  models  in  two  steps;  (i)  the  construction  of 
constitutive  models  which  quantify  the  inherent  hysteresis  and  nonlinearities  in  the  magnetostrictive  material,  and 
(ii)  the  incorporation  of  the  constitutive  relations  into  transducer  models  which  specify  the  device  displacement 
outputs  due  to  specified  current  inputs. 

2.1  Nonlinear  Hysteretic  Constitutive  Models 

we  describe  here  a  domain  wall  model  quantifying  the  hysteretic  relationship  between  an  applied  magnetic  field  H  and 
the  magnetization  M  in  a  Terfenol-D  rod.  The  premise  of  this  model  is  that  one  source  of  hysteresis  in  ferromagnetic 
materials  can  be  attributed  to  the  impediment  of  motion  of  domain  walls  that  are  pinned  at  defects  or  inclusions 
in  the  material.  We  provide  here  those  details  of  the  theory  pertinent  to  the  subsequent  development  of  nonlinear 
adaptive  estimation  techniques  and  refer  the  reader  to  [14,  15]  for  details  regarding  the  development  of  the  theory 
for  general  ferromagnetic  materials  and  [4,  7,  8]  for  the  extension  of  the  theory  to  magnetostrictive  transducers. 

The  domain  wall  model  for  the  magnetization  M  produced  by  an  applied  magnetic  field  H  will  be  developed  in 
three  steps;  (i)  quantification  of  the  anhysteretic  magnetization  Man ,  (ii)  quantification  of  the  irreversible  magneti¬ 
zation  Mirr,  and  (iii)  quantification  of  the  total  magnetization  M . 

Physically,  Man  can  be  interpreted  as  the  equilibrium  magnetization  obtained  in  the  absence  of  material  inclusions. 
Mathematically,  the  anhysteretic  magnetization  can  be  modeled  by  quantifying  the  probability  of  achieving  specified 
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energy  states  through  Boltzmann  relations.  As  detailed  in  [27,  31]  —  or  [28]  for  the  analogous  ferroelectric  model  — 
the  assumption  that  magnetic  moments  can  orient  only  in  the  direction  of  the  applied  field,  or  diametrically  opposite 
to  it,  yields  the  Ising  relation 

Man  =  Ms  tanh  ■  (2) 

Alternatively,  if  one  assumes  that  moments  have  equal  probability  of  orienting  in  any  direction,  one  obtains  the 
Langevin  model 

Man  =  Ms  ^coth  (3) 

which  agrees  with  (2)  through  first-order  terms  —  see  also  [5].  In  both  models,  Ms  denotes  the  saturation  magneti¬ 
zation  and 

He  =  H  +  aM 


quantifies  the  effective  fields  where  a  quantifies  denotes  a  coupling  coefficient, 
relation 


a(T)  = 


HhT 

Tc 


A  Boltzmann  derivation  yields  the 


(4) 


where  Hh  is  a  biasing  field  and  Tc  denotes  the  Curie  temperature.  Whereas  (4)  quantifies  the  temperature-dependence 
of  a  in  certain  operating  regimes,  the  predictions  that  it  provides  should  typically  be  considered  as  qualitative  rather 
than  quantitative  in  nature.  Hence,  the  relation  (4)  is  best  interpreted  as  indicating  that  temperature-dependence 
should  be  accommodated  in  some  manner  —  e.g.,  adaptive  estimation  techniques  —  rather  than  providing  an  accurate 
model  quantifying  the  temperature  behavior. 

The  anhysteretic  relations  can  be  utilized  to  approximate  the  magnetization  at  low  drive  levels  but  they  do  not 
incorporate  the  energy  loss  due  to  the  movement  of  domain  walls  and  thus  can  not  model  high  drive  levels.  To 
quantify  the  total  magnetization,  it  is  necessary  to  incorporate  the  reversible  magnetization  Mrev  and  irreversible 
magnetization  Mirr.  Reversible  magnetization  occurs  when  applied  field  levels  are  sufficiently  small  that  domain 
walls  remain  pinned  while  the  irreversible  component  represents  the  movement  of  pinning  sites  which  occurs  at 
higher  drive  levels.  As  detailed  in  [8,  15],  Mirr  can  be  quantified  as  the  anhysteretic  component  minus  the  losses 
required  to  break  pinning  sites 


-  k 


dMirr 

dHe 


where  k  is  a  measure  of  the  average  energy  required  to  translate  the  domain  wall.  This  relation  can  be  reformulated 
in  terms  of  the  applied  field  as 


dMjrr  _  ~  Man  -  Mirr 
dH  k6  -  a(Man  -  Mirr) 


where  5  =  sign  (dH)  ensures  that  pinning  opposes  change  in  magnetization.  The  switch  <5  is  0  if  dH  >  0  and  M  >  Man 
or  dH  <  0  and  M  <  Man  and  1  otherwise.  This  is  necessary  to  model  the  physical  observation  that  after  a  field 
reversal,  the  changes  in  magnetization  are  purely  reversible  until  the  anhysteretic  value  is  reached. 

The  reversible  magnetization  is  given  by  the  algebraic  relationship 


Mrev  —  c(Man  -\fjrr ) 


where  c  is  a  material  parameter  which  quantifies  the  reversibility  of  the  material.  The  total  magnetism  is  the  sum  of 
the  reversible  and  irreversible  magnetizations  which  yields 


M  =  (1  —  c)Mirr  +  cMan. 


The  relation  (5)  can  be  reformulated  as  a  differential  equation 


dM 

~dH 


F(H ,  M) 


where 


M(H0)  =  M0 


F(H,  M) 


1 


1  +  ca 


dMarl 

dH 


'?  Man-M  dMa, 

kS  —  a(Man  —  M)  dH 


(5) 

(6) 
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with  a  =  .  The  anhysteretic  magnetization  Man  can  be  taken  as  either  the  Ising  model  (2)  or  the  Langevin 

model  (3). 

The  anhysteretic  behavior  modeled  by  the  Ising  relation  (2)  and  total  magnetization  specified  by  (5)  or  (6) 
are  illustrated  in  Figure  2.  The  parameter  values  are  a  =  7012  A/m,  k  =  400  A/m,  a  =  —0.1,  Ms  =  7.65  x 
105  A/m  and  c  =  0.18.  which  were  determined  in  [8]  through  a  least  squares  fit  to  data  measured  under  isothermal 
conditions  at  room  temperature.  We  point  out  that  the  model  exhibits  a  nonlinear  dependence  on  the  parameters 
k  and  c  arising  in  Mjrr  and  Mrev  as  well  as  the  parameters  Ms,  a  and  a  arising  in  Man.  This  necessitates  the 
development  of  nonlinear  adaptive  identification  techniques  to  update  parameter  values  online  to  accommodate 
changing  environmental  conditions  such  as  temperature. 

2.2  Transducer  Models 

The  models  (5)  or  (6)  quantify  the  nonlinear  hysteretic  relationship  between  the  imposed  field  and  the  resulting 
magnetization  in  a  magnetostrictive  material.  However,  the  magnetoelastic  coupling  which  provides  the  materials 
with  magnetostrictive  properties  has  not  yet  been  addressed.  In  applications  where  stress  anisotropies  dominate 
crystalline  anisotropies,  experimental  data  [8]  suggests  a  quadratic  dependence  of  strain  on  magnetization.  This 
can  be  theoretically  justified  by  noting  that  in  present  transducer  designs,  the  prestress  levels  are  sufficient  to  allow 
stress  anisotropies  to  dominate  crystalline  anisotropies  so  that  the  strains  are  primarily  due  to  the  quadratic  rotation 
processes.  To  incorporate  this  quadratic  dependence  the  stress-induced  in  the  Terfenol-D  rod  is  specified  as 

a  =  Ye-  YjM2  (7) 

where  Y  denotes  Young’s  modulus  and  7  is  a  magnetoelastic  coupling  coefficient.  This  relationship  quantifies  the 
linear  Hooke’s  relationship  between  stresses  a  and  strains  e  in  the  Terfenol-D  rod  as  well  as  the  nonlinear  hysteretic 
dependence  of  the  stresses  on  applied  fields  H  through  the  magnetization  M. 

The  relation  (7)  quantifies  the  local  elastic  behavior  of  the  Terfenol-D  rod  but  does  not  include  spatial  dependence 
nor  internal  damping  in  the  transducer.  We  now  develop  models  which  specify  the  displacements  and  forces  generated 
by  the  magnetostrictive  transducer  due  to  an  input  current  I(t).  A  partial  differential  equation  model  developed  in 
[8]  and  the  corresponding  ordinary  differential  equation  (ODE)  system  derived  through  a  finite  element  discretization 
are  first  summarized.  We  then  demonstrate  that  the  ODE  system  can  be  adequately  approximated  by  a  second-order 
scalar  differential  equation  for  transducers  which  employ  flux  shaping  techniques  to  minimize  end  effects. 

The  left  end  of  the  Terfenol-D  rod  (x  =  0)  is  assumed  to  be  fixed  while  the  right  end  (x  =  L)  is  constrained  by 
a  damped  oscillator  and  has  an  attached  point  mass,  as  depicted  in  Figure  3.  The  Kelvin- Voigt  damping  coefficient 
and  density  of  the  Terfenol-D  rod  are  denoted  by  cd  and  p,  and  the  point  mass  is  denoted  by  Ml.  The  damping 
spring  has  stiffness  A7  and  damping  coefficient  cl ■  Strains  in  the  rod  are  specified  in  terms  of  the  longitudinal  rod 
displacement  w  by  e  =  ^ . 

The  direct  use  of  the  constitutive  relation  (7)  will  yield  an  undamped  model  for  the  Terfenol-D  rod.  We  assume 
that  stress  is  proportional  to  a  linear  combination  of  strain  and  strain  rate  to  incorporate  Kelvin- Voigt  damping  . 


Applied  Field 

(a) 


(b) 


Figure  2:  (a)  Anhysteretic  magnetization  given  by  (2),  and  (b)  total  magnetization  specified  by  (5)  or  (6). 
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+u 


x=L 


Figure  3:  Rod  model  for  Terfenol-D  transducer. 


Thus,  the  stress  at  any  point  x,  0  <  x  <  L,  is  given  by 


,  .  „rdw,  .  d2w  ,  .  ,,  . 

a(t,x)  =  Y—(t,x)  +  cD-^^(t,x)  —  Y^M2(t,x). 


(8) 


The  relation  (8)  is  identical  to  the  relation  obtained  in  [8]  if  the  coupling  coefficient  is  defined  as  7  =  A S/Ms 
where  Xs  and  Ms  respectively  denote  the  saturation  magnetostriction  and  saturation  magnetization,  respectively. 
Balancing  forces  yields 

^  d2w  _  dNtot  ,  , 

PA  dt 2  dx  (  ^ 

where  A  is  the  cross  sectional  area  of  the  Terfenol-D  rod  and  the  force  resultant  is  specified  by 

dw  d2  xu 

Ntot(t,x)  =  YA—(t,x)  +  cDA—-—-(t,x)  -  Y A^M(t,x)2 . 

OXJJ  oxot 

To  obtain  appropriate  boundary  conditions,  we  first  note  that  w(t,  0)  =  0.  Balancing  forces  at  x  =  L  yields 

dw  d2xu 

Ntot(t,L)  =  —kLw(t,  L)  -  cL  —  (t,L)  -  MLQ-^(t,L). 
dxu 

Initial  conditions  are  taken  to  be  iu(0,a;)  =  0  and  — —  (0,  a;)  =  0. 

ox 

As  detailed  in  [8],  formulation  of  the  model  (9)  in  weak  form  and  spatial  discretization  utilizing  linear  finite 
elements  yields  the  vector  ODE  system 


MS  +  C£+K.x  =  -iB  [. M2(H)}  (t) 
a?(0)  =  So,  S(  0)  =  S  i 


(10) 


where  x(t)  €  1ZN  specifies  the  states,  M  €  lZNxN ,  C  £  TZNxN ,  K,  £  7 ZNxN  respectively  denote  the  mass,  damping 
and  stiffness  matrices,  and  B  £  7 ZN  contains  integrated  basis  functions  resulting  from  the  inputs.  For  general  systems, 
large  N  (e.g.,  N  =  32)  may  be  required  to  achieve  convergence.  This  may  limit  the  speed  at  which  model-based 
control  designs  can  be  implemented. 

For  transducer  design  in  which  the  permanent  magnet  is  constructed  to  minimize  end  effects,  measurements  with 
a  Hall  probe  indicate  that  nearly  uniform  fields 


H(t,x)  =  nl(t )  (11) 

are  achieved  along  the  length  of  the  rod  here  n  denotes  the  number  of  coils  per  unit  length  in  the  solenoid.  In 
such  cases,  each  finite  element  section  of  the  rod  reacts  identically  to  the  uniform  magnetic  field  which  motivates  the 
consideration  of  a  lumped  scalar  model 


mx  +  cx  +  kx  =  £  \M2(H)\  (t) 
i(0)  =  xq  x(0)  =  So 

where  m,  c,  k  and  £  are  parameters  which  must  be  estimated  for  a  given  transducer  construction. 


(12) 
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Figure  4:  Comparison  of  the  lumped  model  (13)  with  the  finite  element  model  (10)  linearized  about  the  bias  level 
( H0,M0 )  with  (a)  N  =  2,  and  (b)  N  =  256. 


Experiments  also  indicate  that  for  moderate  drive  levels  about  a  bias  ( Hq ,  Mo),  stresses  exhibit  an  approximately 
linear  dependence  on  the  magnetization  but  a  nonlinear  and  hysteretic  dependence  on  the  field  H  or  current  I  to  the 
solenoid.  This  observation  motivates  us  to  linearize  about  the  bias  level  Ma  to  obtain  the  final  model 

x  +  cx+  kx  =  £[M(H)]  (t) 
a;(0)  =  Xq  i(0)  =  xq- 

The  scalars  c,  k  and  £  are  determined  through  a  fit  to  the  Galerkin  approximation  (10)  or  data  from  the  physical 
device.  A  least  squares  fit  to  the  Terfenol-D  data  reported  in  [8]  yielded  the  parameter  values  c  =  7.8899  x  103, 
k  =  6.4251  x  107  and  £  =  1.3724  x  10-2.  The  accuracy  of  the  lumped  model  (13)  for  transducer  designs  having 
minimal  end  effects  is  illustrated  in  Figure  4  where  it  is  compared  with  the  finite  element  model  (10)  linearized  about 
the  bias  level  (i7o,Mo).  The  accuracy  of  the  finite  element  model  for  characterizing  the  physical  device  is  reported 
in  [7,  8].  Finally,  we  note  that  while  the  model  (13)  exhibits  a  linear  dependence  on  M,  it  retains  the  fully  hysteretic 
and  nonlinear  dependence  on  H  and  I  through  the  relations  (6)  and  (11). 

3  Nonlinear  Adaptive  Parameter  Estimation  Algorithms 

The  previous  section  summarized  the  PDE  model  developed  in  [8]  and  discussed  the  development  of  a  second  order 
scalar  ODE  model  (13)  quantifying  the  tip  displacement  generated  by  an  applied  field  or  current  in  the  Terfenol-D 
transducer  depicted  in  Figure  la.  The  material  parameters  in  the  magnetization  model  (6)  were  initially  identified 
under  the  assumption  of  a  constant  operating  temperature.  However,  in  industrial  applications  the  assumption  of 
fixed  temperature  is  often  invalid  due  to  Joule  heating  resulting  from  eddy  currents  and  ohmic  heating  in  the  solenoid. 
This  necessitates  the  development  of  nonlinear  adaptive  identification  techniques  to  update  slowly  varying  material 
parameters  to  accommodate  the  effects  of  changing  temperatures. 

In  particular,  the  model  parameter  a  in  the  nonlinear  hysteretic  magnetization  model  (6)  has  been  observed  to 
reflect  changes  in  behavior  due  to  temperature  as  noted  in  (4)  and  the  subsequent  discussion.  Figure  5  illustrates 
the  effect  of  temperature  changes  on  the  H-M  relationship  quantified  by  the  model  (6)  for  two  different  values  of  the 
parameter  a.  In  this  initial  development,  we  focus  on  the  development  of  nonlinear  adaptive  parameter  estimation 
techniques  to  update  a  to  accommodate  temperature-dependencies. 

To  accomplish  this,  we  consider  the  theory  in  [17,  18]  and  develop  modifications  required  for  the  magnetostrictive 
transducer  model  employed  here.  To  begin,  we  will  summarize  the  relevant  aspect  of  the  nonlinear  adaptive  parameter 
estimation  algorithms  discussed  in  [17,  18].  Since  the  algorithms  developed  in  [17,  18]  identify  nonlinearly  occurring 
parameters  in  first-order  differential  equations,  we  must  subsequently  extend  these  concepts  in  order  to  develop  a 
nonlinear  adaptive  parameter  estimation  algorithm  appropriate  for  second-order  ODE  systems  such  as  the  model 
(13)  developed  to  quantify  transducer  behavior  in  response  to  nonlinear  and  hysteretic  inputs. 
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Figure  5:  The  H-M  relation  quantified  by  (6)  for  two  different  values  of  a: 


(-)  a  =  7012,  (-  -)  a  =  6300. 


3.1  Scalar  Case 

The  nonlinear  parameterization  techniques  developed  in  [17,  18]  are  based  on  the  assumption  that  all  states  are 
available  in  systems  of  the  form 

V  =  -A  y  +  vf(u(t),6) 

where  A  >  0  and  v  are  scalars,  9  £  0  is  an  unknown  parameter  and  0  is  the  bounded  region  in  Rm.  The  function  / 
is  taken  to  be  a  scalar- valued  nonlinear  function  of  the  input  u(t).  For  example,  in  the  numerical  results  presented 
in  Section  4,  /  is  taken  to  be  the  magnetization  model  (6),  u(t)  is  the  applied  field  or  current,  and  6  is  taken  to  be 
the  nonlinearly  occurring  parameter  a.  As  motivated  by  [17,  18],  we  consider  the  estimation  algorithm 


V 

V 


Ve 

e 


-A 


y  —  e  sat 


+  f(u,  9)  —  Csat 


y  -  y 


y  —  e  sat  [  - 


where  e  >  0,  sat(-)  is  a  saturation  function  defined  by 


sat(x) 


1, 
x , 
-1, 


X  >  1 

1*1  < 1 

X  <  —  1 


and  u*  and  <j>*  are  the  solution  of 


v*  =  min  max  J(9,d>) 

06 1"  060 

d>*  =  arg  min  max  J(9,d>) 

0eRm  060 


Wt) 


/(M)-/(M)-0T(0-0) 


(14) 


(15) 


We  note  that  when  |y|  <  e,  the  adaptation  of  the  parameters  stops,  imposing  what  is  termed  a  dead-zone.  The 
method  will  continue  to  adapt  the  parameters  until  the  magnitude  of  the  error  y  is  less  than  the  given  e. 

We  consider  the  min-max  algorithm  (15)  to  handle  the  regions  of  nonconvexity  of  /  where  the  gradient  method 
is  insufficient  as  discussed  in  Section  1.  The  use  of  a  tuning  error  ye  rather  than  a  tracking  error  y  ensures  continuity 
of  the  adaptation  as  does  the  use  of  a  saturation  function  rather  than  a  signum  function  [18]. 

Defining  9=9  —  9  and  x  =  \y1  9T]T,  it  can  be  shown  that  the  system  (14)  is  stable  about  y  =  0  by  proving  that 

V  =  y2e  +  92 
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is  a  Lyapunov  function.  To  illustrate  the  manner  through  which  this  is  established,  we  summarize  details  from  the 
original  proof  in  [18]. 

Summary  of  Proof  from  [18]:  We  first  note  that 


1/  =  2  yeye  +  289. 

If  \v\  <  e,  then  ye  =  0  so  it  follow  that  V  =  0.  We  then  need  to  show  that  V  <  0  if  \y\  >  e.  We  can  represent  V  as 


y  =  2  ye 


-A  [y  -  e  sat  -  +  f(u,  9)  -  i/*sat  -  +  ay  -  f(u ,  9) 


-  2 9ye<j>* 


=  -2Aye  ye  +  2 y, 
=  — 2A  y\  +  2  y, 


f(u ,  0)  —  f(u ,  0)  —  0(f)*  —  v*sat 
f(u,9 )  —  f(u,w)  —  9(j)*^j  —  i/*  sat 


If  y  >  0,  then  sat  =  1  so  we  must  have 


v  >  sat  (  - 


f(u,9)  -  f(u,8)  -  wcj>* 


for  all  9  £  0. 


This  implies  that  we  can  let 


*  ,  ,  V 

v  =  max  sat  |  — 
060 


f(u,  9)  -  f(u,  8)  -  wcf*  for  any  cf* 


The  inequality  (16)  is  satisfied  by  the  definition  of  <jf  and  v*  so  that  V  <  0. 
If  y  <  0,  then  sat  =  —1  so  we  must  have 


v*  >  - 


or 


v  >  sat 


f(u,  9)  —  f(u,  9)  —  6( 
f{u,9)  -  f{u,0)  -  w<j) 


for  all  9  £  0 


for  all  9  £  0. 


Again  we  can  let 


*  4.  i  y 

v  =  max  sat  — 

060 


f{u,8)  -  f(u,9)  -  w(j)* 


for  any  cf* . 


By  the  definition  of  (f> *  and  v*  the  inequality  (17)  is  satisfied,  hence  V  <  0. 


(16) 


(17) 


□ 


To  implement  the  method  proposed  in  the  system  (14),  it  is  necessary  to  solve  the  min-max  problem  (15).  To  do 
so,  we  must  construct  a  concave  cover  F{8)  and  convex  cover  F_(&)  where  the  covers  satisfy 

p(o) >f-7  m <f-T 

with  /  =  f(u,9).  The  following  definitions  and  construction  are  summarized  from  [18]. 

Definition  1:  A  point  9°  £  9C  if  9°  £  0  and 

V/flo(0-0°)>/-/0  (18) 

df 


where  V/go  = 


89 


and  f°  =  f(<j),  9°). 


00 


Definition  2:  9C  =  9C  n  0  where  9C  is  the  complement  of  9C. 
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Figure  6:  Convex  cover  F_(0)  and  concave  cover  F{9)  of  f(u,9). 


If  /  is  not  concave  on  0,  then  0C  is  given  by  9C  =  {912,934, . . . ,  0mn}  where  9^  =  [9l ,  0J],  >  9l ,  are  the  regions 

where  /  is  not  concave.  Using  Definitions  1  and  2,  the  concave  cover  of  /  —  /  on  0  can  be  subsequently  constructed 
as  ^ 

_  f  /  —  /,  for  all  9  £  9C 

^)=  ..  ..  ..  -  (19) 

I  uj1-W  +  clP  for  all  9  £  913  £  9C 


where 


.‘J  = 


P  -  P 


cij  =  P  -  f  -  Pj9\  f  = 


Qi  -  9i 

Similarly,  we  construct  a  convex  cover  of  /  —  /  by  defining 

9V  =  {0°  |  Vfgo{9-9°) 

L  =  in© 


E(0)  = 


f-f, 

uij9  +  cij, 


for  all  9  £  9V 
for  all  9  £  9 v  £  9V  . 


(20) 


An  illustrative  example  of  a  convex  and  concave  cover  of  a  function  /  is  depicted  in  Figure  6. 

Once  we  have  constructed  F_{6)  and  F{9):  a  closed  form  solution  to  the  min-max  problem  (15)  is  given  by 


m 

I  Vfs,  if  9  £9c 
\  Pj,  if  9  £  9ij  £  9C 

- m 

V/g,  if  9  £  9V 
u>ij,  if  9  £  9l]  £  L 


if  Ve  >  0 


(21) 


if  S/e  <  0  . 


A  proof  that  (21)  is  the  solution  to  (15)  can  be  found  in  [18]. 

Having  established  the  stability  of  the  adaptation  method  by  the  Lyapunov  function,  we  now  seek  sufficient  condi¬ 
tions  which  guarantee  uniform  asymptotic  stability  of  the  system  (14).  We  summarize  the  condition  for  convergence 
of  the  parameters  and  tracking  error  in  the  following  theorem  presented  from  [17]. 

Theorem  1:  If  for  every  t\  >  to,  there  exists  To,  eo,  So,  and  a  subinterval  fe,  £2  +  <$o]  €  [£i,  £1  +  To]  such  that 

pt2+So 


P 


ft  2 


f(u,  9(tf))  —  f(u,  9)  dr  >  2e  +  e0|  |6»(£2)||  , 


(22) 


then  the  origin  y  =  0  is  uniformly  asymptotically  stable. 
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In  Theorem  1,  (3  =  1  if  f(u,9)  is  convex  and  (3  =  —  1  if  /(it,  0)  is  concave.  We  note  several  differences  between 
this  condition  and  the  condition  for  a  linear  parameterization.  First,  the  sign  of  the  integral  is  important.  The  sign 
is  not  strictly  determined  by  f(u,9)  —  f{u,9)  but  also  by  the  convexity  or  concavity  of  /  as  indicated  by  f3.  This 
coupling  arises  from  the  min-max  algorithm  and  is  necessary  but  not  sufficient  to  ensure  that  the  method  will  leave 
the  dead  zone  |j/|  <  e.  The  integral  must  also  be  sufficiently  large  to  leave  the  deadzone,  which  necessitates  the  term 
incorporating  e  on  the  right  hand  side  of  (22). 

Whereas  Theorem  1  places  the  excitation  conditions  on  /,  we  wish  to  derive  conditions  on  u(t)  since  we  have 
some  freedom  when  choosing  u(t).  Theorem  1  does  not  give  conditions  on  the  input  u  to  satisfy  the  inequality  (22) 
nor  does  it  guarantee  that  such  an  input  exists.  Inequality  (22)  includes  two  components.  First,  the  magnitude  of  the 
integrand  must  be  sufficiently  large.  For  a  large  parameter  error,  the  input  must  be  such  that  the  difference  between 
the  function  evaluated  at  the  actual  parameter  and  the  parameter  estimate  is  adequately  large;  therefore  we  desire 
an  input  signal  which  drives  the  function  /  to  a  level  where  a  change  in  the  parameter  is  most  noticeable.  Secondly, 
the  integral  must  be  the  same  sign  as  f3.  This  condition  states  that  if  /  is  convex,  then  the  integrand  should  be 
positive  and,  conversely,  if  /  is  concave,  then  the  integrand  should  be  negative.  The  min-max  feature  of  the  algorithm 
provides  stability,  but  an  acceptable  input  must  be  used  to  guarantee  parameter  convergence.  Parameter  convergence 
is  ensured  by  updating  using  the  gradient  information  and  we  must  pick  an  input  signal  accordingly. 

To  ensure  parameter  convergence,  we  summarize  the  persistent  excitation  conditions  developed  in  [17]  for  the 
input  it. 

Persistent  Excitation  Conditions: 

3.1(a).  For  a  given  9 ,  it  must  reverse  the  sign  of  the  integrand  of  (22)  while  keeping  the  convexity /concavity  of  / 
fixed, 
or 

3.1(b).  For  a  given  (9,  u  must  reverse  the  convexity /concavity  of  /,  while  preserving  the  sign  of  the  integrand  in 
equation  (22). 

A  numerical  example  illustrating  these  conditions  is  provided  in  Section  4. 

3.2  Vector  System 

The  nonlinear  identification  algorithm  developed  in  [17,  18]  and  summarized  in  Section  3.1  is  for  scalar  differential 
equations  and  hence  is  not  directly  applicable  to  higher-order  differential  equations  or  systems  of  differential  equations 
of  the  type  arising  in  physical  applications  such  as  the  quantification  of  hysteresis  in  a  magnetostrictive  transducer. 
In  this  section,  we  extend  the  nonlinear  parameter  identification  algorithm  to  systems  of  differential  equations.  We 
develop  the  analysis  in  the  context  of  the  motivating  second-order  system  but  note  that  it  is  sufficiently  general  to 
include  higher-order  systems. 

To  extend  the  method  to  vector  systems,  we  must  redefine  several  of  the  variables  from  the  scalar  case.  Since 
we  want  to  use  the  solution  to  the  min-max  problem  (15),  we  must  ensure  that  we  do  not  affect  that  aspect  of  the 
formulation. 

Consider  the  identification  of  nonlinearly  occurring  parameters  in  the  matrix  system 

y  =  Ay  +  Bf(u,9)  (23) 

where  A  £  R2x2  and  B  £  R2xl.  We  consider  the  adaptation  algorithm 

/(it,  9)  —  a*sat 

y  =  y-y  (24) 

Q  =  ~iyt<t>* 

ye  =  B2p2y  l  +  B2p3y2 


V  =  Ay  +  B 


11 


where  P2  and  p$  are  scalars  to  be  determined,  7  >  0  is  a  scalar  gain,  B2  is  the  (2,1)  component  of  vector  B,  and  a 
and  (f>*  are  solutions  of  the  min-max  problem 


a*  =  min  max  q(9 ,  <b) 
0eRm  see  v 


d>*  =  arg  min  max  q(9 ,  d>) 

4>mm  Bee  y 


(25) 


9(8  A) 


f{u,  8)  —  f{u,  9)  —  <j>T(d  —  9) 


It  is  important  to  note  that  the  solution  to  the  min-max  problem  (25)  is  the  solution  to  (15)  so  this  aspect  of  the 
adaptive  estimation  is  the  same  as  the  scalar  case.  A  notable  differences  between  the  matrix  system  method  (25)  and 
the  scalar  method  (15)  is  the  definition  of  yt.  The  alteration  of  ye  is  necessary  for  the  Lyapunov  candidate  utilized 
for  the  vector  system. 

The  following  theorem  ensures  the  globally  asymptotic  stability  of  the  adaptive  parameter  estimation  algo¬ 
rithm  (24). 

Theorem  2:  If  the  vector  differential  equation  (23)  is  globally  asymptotically  stable,  the  adaptive  parameter  estima¬ 
tion  algorithm  (24)  is  globally  asymptotically  stable. 

Proof:  We  consider  the  Lyapunov  function 

V  =  'yyTPy  +  02  (26) 

where  P  is  a  symmetric  positive  definite  matrix.  Taking  the  derivative  of  (26)  yields 


V  =  "fyTPy  +  "fy1  Py  +  2  69 


with 


y  =  Ay  +  B 


f  -  f-  sat  (y) 


Employing  the  definitions  of  y  and  9 ,  we  can  write 


V  =  7  (i  fATPy  +  fPAy)  +  27  ( 


/  -  /  -  a*  sat  (  —  ; 


(?) 


BTPy  -  ije 


')■ 


(27) 


Since  A  is  globally  asymptotically  stable,  for  any  symmetric  positive  definite  matrix  Q  there  exists  a  symmetric 
positive  definite  matrix  P  such  that  AT P  +  PA  =  —Q  (e.g.,  see  [2]).  Denoting 


Pi  P2 
P2  P3 


we  have  BT Py  =  u>p2yi  +  u>p^y2  =  ye ■  Substituting  this  into  (27)  yields 

V  =  - 7  ( yTQy )  +  27 ye  (f  -  f  -  8<j>*  -  a*sat  (y )  )  ■ 


The  proof  that  V  <  0  is  completed  by  utilizing  the  definitions  of  a*  and  <f>*  as  the  solutions  of  (25)  in  a  manner 
similar  to  that  employed  in  the  proof  in  the  scalar  case.  Note  that  the  parameters  P2  and  p%  can  be  determined  for 
any  symmetric  positive  definite  choice  of  Q.  □ 


Notable  differences  between  the  vector  parameter  estimation  algorithm  and  the  scalar  algorithm  include  the  lack 
of  a  dead-zone  and  the  omission  of  the  esat(|)  term  in  the  definition  of  y.  These  alterations  are  due  to  the  form  of 
the  Lyapunov  function. 


4  Numerical  Examples 

Scalar  and  vector  examples  are  provided  to  demonstrate  properties  of  the  nonlinear  adaptive  parameter  estimation 
algorithms.  In  both  examples,  the  nonlinear  function  is  taken  to  be  the  (6)  for  the  magnetization  M  with  the  input  u 
being  the  applied  field  H .  As  stated  previously,  the  nonlinearly  occurring  parameter  a(T)  in  the  magnetization  model 
reflects  changes  in  the  transducer  behavior  due  to  temperature  changes.  We  make  the  assumption  that  the  value 
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Figure  7:  Control  trajectory  for  the  milling  process  and  regime  during  which  adaptive  estimation  occurs. 

of  a  is  slowly  varying  so  that  the  value  can  assumed  to  be  fixed  during  the  nonlinear  parameterization  estimation 
process.  During  operation,  the  temperature  of  the  transducer  may  change  and  the  corresponding  value  of  a  can  be 
re-identified  during  the  operation  of  the  transducer. 

To  provide  a  specific  operating  regime,  we  consider  the  use  of  the  magnetostrictive  transducer  depicted  in  Figure  la 
for  milling  out-of-round  automotive  products  at  speeds  of  3000  rpm  and  tolerances  of  1  //m  as  described  in  [26] .  The 
specified  trajectory  for  the  cutting  head  attached  to  the  Terfenol-D  rod  is  depicted  in  Figure  7.  The  cutting  head 
starts  from  rest,  is  moved  to  a  prescribed  position  adjacent  to  the  ingot,  and  then  is  cycled  periodically  to  produce 
the  desired  shape  in  the  automotive  product.  It  is  during  the  initial  positioning  stage  that  the  nonlinear  adaptive 
identification  algorithms  are  utilized  to  accommodate  thermally-induced  changes  in  the  transducer  behavior,  and  it 
is  only  during  the  subsequent  periodic  cutting  stage  that  the  accuracy  of  1  (im  must  be  maintained. 

For  both  the  scalar  and  vector  examples,  the  parameters  in  the  magnetization  model  (6)  are  taken  to  be  k  = 
400  A/m,  a  =  —0.1,  Ms  =  7.65  x  105  A/m,  c  =  0.18  and  the  true  value  of  a  is  taken  as  7012  A/m.  Also,  both 
examples  assume  that  a  is  bounded  with  a  £  [6300,  7300]  and  a(0)  =  6990.  The  bounds  imposed  on  a  can  be 
determined  from  experiments  under  various  temperature  regimes. 

4.1  Scalar  Case 

We  illustrate  first  the  scalar  case.  The  dynamics  of  the  system  are  specified  by  the  differential  equation 

V  =  -A  y  +  M(H,a) 

with  the  scalar  A  is  taken  to  be  100  and  M  specified  as  the  solution  to  (5)  or  (6). 

One  issue  in  the  adaptive  parameter  estimation  algorithm  is  the  construction  of  inputs  u(t)  =  H(t )  which  will 
guarantee  persistent  excitation.  The  conditions  imposed  for  excitation  prompt  the  use  of  a  signal  that  does  not 
cause  the  function  M  to  change  signs.  It  is  observed  that  a  monotonically  increasing  input  provides  persistent 
excitation  as  well  as  quick  convergence.  The  input  field  H(t)  was  chosen  to  be  a  linear  function  which  drives  the 
magnetization  to  levels  where  variations  in  the  parameter  a  are  most  significant  as  well  as  realistically  simulating  the 


(a)  (b) 

Figure  8:  (a)  Value  of  the  integrand  in  (22),  and  (b)  resulting  value  of  M(u,a)  provided  by  (5)  or  (6). 
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(a)  (b) 

Figure  9:  (a)  Parameter  estimate  a  and  true  value  a  =  7012  A/m,  and  (b)  tracking  error  y  for  scalar  case. 


initialization  procedure  for  the  milling  application.  Figure  8a  illustrates  the  integrand  of  (22)  for  a  given  value  of  6 
to  demonstrate  that  Persistent  Excitation  Condition  3.1(b)  is  met.  The  integrand  remains  positive  while  switching 
the  convexity /concavity  of  the  function  M  as  seen  in  Figure  8b. 

Figure  9  illustrates  the  capability  of  the  scalar  nonlinear  parameter  estimation  method  to  accurately  estimate  the 
unknown  parameter  a.  Figure  9a  depicts  the  evolution  of  the  parameter  estimate  a  which  quickly  converges  from  the 
intial  value  of  a  =  6990  A/m  to  the  true  value  of  a  =  7012  A/m.  The  tracking  error  y  is  plotted  in  Figure  9b.  The 
speed  of  convergence  of  the  estimation  algorithm  is  important  since  this  identification  method  could  potentially  be 
combined  with  a  control  technique  for  online  estimation  prior  to  control.  We  note  that  the  time  scales  in  Figures  8 
and  9  are  scaled  to  that  in  Figure  7  -  the  convergence  rates  of  both  the  integrand  of  (22)  and  the  estimation  algorithm 
are  determined  by  the  time  needed  to  cycle  the  magnetization  which  can  be  modified  according  to  the  application. 


4.2  Vector  System 


We  illustrate  here  the  nonlinear  adaptive  parameter  algorithm  for  a  vector  system  as  developed  in  Section  3.2.  The 
system  we  consider  is  the  damped  spring  mass  system  (13)  used  to  model  the  transducer  dynamics  in  response 
to  nonlinear  and  hysteretic  inputs  with  the  coefficients  specified  in  Section  2.2.  The  values  of  the  domain  wall 
model  parameters  are  those  specified  in  Section  4.1.  Again,  the  parameter  a  is  to  be  determined  because  of  its 
temperature-dependence.  The  gain  7  is  taken  to  be  5  x  IO10  and  P  is  found  utilizing  the  matrix 


Q 


1  X  IO10  0 
0  lx  IO"8 


We  employ  the  same  applied  magnetic  field  as  in  the  scalar  example. 

The  convergence  of  the  estimate  A  for  the  vector  system  is  illustrated  in  Figure  10  while  Figure  11  illustrates  the 
tracking  error  of  the  adaptive  system  for  the  vector  system  in  which  y  =  [i/1,2/2]  denotes  the  tracking  error  in  the 
states  y  =  [221,2/2]- 


Figure  10:  Parameter  estimate  a  and  true  value  a  =  7012  A/m  for  the  vector  system. 
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(a) 


(b) 


Figure  11:  Tracking  errors  (a)  in  the  displacement  jq,  and  (b)  velocity  i/2  for  the  vector  system. 


There  are  a  number  of  numerical  issues  associated  with  the  implementation  of  the  estimation  algorithms  for  the 
scalar  or  vector  systems.  Inaccuracy  in  the  solution  can  cause  the  value  of  y  to  have  discontinuous  jumps  from 
positive  to  negative  values.  This  phenomena  causes  the  min-max  solution  to  jump  between  utilizing  the  convex  cover 
and  concave  cover  which,  in  turn,  can  cause  highly  oscillatory  behavior  in  the  parameter  updates.  The  models  must 
be  solved  numerically  and  since  implicit  methods  are  not  feasible  because  of  the  unknown  forcing  function  at  the 
next  time  step,  the  time  step  must  be  sufficiently  small  to  ensure  accurate  solutions  of  the  model.  These  issues  have 
been  addressed  numerically  and  the  experimental  implementation  of  the  algorithm  is  under  present  investigation. 


5  Concluding  Remarks 

The  nonlinear  parameter  identification  technique  developed  in  [17,  18]  has  been  shown  to  be  effective  for  identifying 
and  updating  a  temperature-dependent  parameter  in  the  nonlinear  domain  wall  model  incorporated  as  the  forcing 
function  in  a  scalar  ODE.  We  have  also  extended  this  parameter  identification  method  to  apply  to  second-order 
systems.  In  this  latter  framework,  this  adaptive  parameter  estimation  algorithm  can  be  used  to  estimate  and 
update  parameters  in  a  wide  range  of  hysteretic  systems  including  SMAs,  ferroelectrics  and  ferromagnetics.  Due 
to  the  persistent  excitation  conditions  placed  on  the  input  to  the  hysteresis  model,  there  are  potential  difficulties 
in  employing  the  scalar  method  as  a  closed-loop  nonlinear  adaptive  control  method  by  considering  the  nonlinear 
function  in  the  dynamics  to  be  the  composition  of  an  approximate  inverse  model  and  the  hysteresis  system  as  was 
done  linearly  in  [33,  34].  The  vector  method  does  not  require  persistent  excitation  conditions  and  hence  it  may  be 
amenable  to  incorporation  in  adaptive  control  methods  utilizing  inverse  models. 

For  control  designs  utilizing  this  parameter  identification  method,  one  can  identify  the  temperature-dependent 
parameter  a  in  the  hysteresis  model  while  the  transducer  is  being  initialized.  Once  the  hysteresis  model  is  completely 
determined,  a  partial  inverse  compensation  method,  e.g.,  [20]  or  [26],  could  be  applied.  The  method  in  this  paper 
can  subsequently  be  employed  to  periodically  re-identify  the  parameters  in  the  hysteresis  model  to  ensure  accurate 
modeling  as  the  parameters  slowly  vary  due  to  the  internal  heat  generated  by  either  the  solenoid  or  surrounding 
media.  This  may  be  more  efficient  and  accurate  than  current  methods  of  parameter  identification,  such  as  least 
square  fits  to  data,  and  can  be  performed  online.  The  next  step  would  be  to  apply  this  method  to  multiple  unknown 
parameters,  i.e. ,  6  G  Rm.  This  method  potentially  could  provide  a  new  method  to  identify  the  full  set  of  parameters, 
{a,  fc,  c,  a,  Ms},  simultaneously.  While  the  method  can  theoretically  accommodate  multiple  parameter  identifica¬ 
tion,  there  are  numerical  difficulties  in  identifying  the  regions  of  concavity  and  convexity.  In  the  one  parameter  case, 
we  need  only  to  identify  linear  functions  to  cover  the  nonconvex  or  nonconcave  sections  of  the  nonlinear  function 
whereas  in  a  multiple  parameter  identification  bounding  hyperplanes  are  required. 
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